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The Geometric Phase Appears in the Ultracold Hydrogen Exchange Reaction 
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Quantum reactive scattering calculations for the hydrogen exchange reaction H + H 2 (u = 4, 
j = 0) —)■ H + H 2 (u/ j') and its isotopic analogues are reported for ultracold collision energies. Due 
to the unique properties associated with ultracold collisions, it is shown that the geometric phase 
effectively controls the reactivity. The rotationally resolved rate coefficients computed with and 
without the geometric phase are shown to differ by up to four orders of magnitude. The effect is 
also significant in the vibrationally resolved and total rate coefficients. The dynamical origin of the 
effect is discussed and the large geometric phase effect reported here might be exploited to control 
the reactivity through the application of external fields or by the selection of a particular nuclear 
spin state. 


In the Born-Oppenheimer description of molecules the 
electronic Schrddinger equation is solved to obtain an ef¬ 
fective potential energy surface (PES) which is then used 
in the solution of the nuclear motion Schrddinger equa¬ 
tion. The electronic PES often becomes degenerate with 
an excited electronic state resulting in a conical intersec¬ 
tion (Cl). As noted long ago by Longuet-Higgins[lj and 
Herzberg and Longuet-HigginSjQ the electronic wave 
functions associated with a Cl change sign for any nu¬ 
clear motion pathway which encircles the Cl (i.e., they 
are double-valued). The electronic sign change implies 
that a corresponding sign change must also occur on 
the nuclear motion wave function. Mead and TruhlarQ 
showed that this can be accomplished by including an 
effective vector potential in the nuclear motion Hamil¬ 
tonian. MeadQ originally referred to this effect as the 
“Molecular Aharonov-Bohm” effect but it is now com¬ 
monly referred to as the “geometric phase” or “Berry’s 
phase” effect. [ 1 ,@ 

The most studied of all chemical reactions is the hydro¬ 
gen exchange reaction H -|- H 2 —>■ H -|- H 2 and its isotopic 
analogues H -|- HD o D -|- H 2 and D -|- HD o H -|- D 2 . 
The H 3 system exhibits a Cl between the ground and first 
excited electronic states for equilateral triangle {i.e., D^y) 
geometries. As first predicted by Mead, the GP alters the 
relative sign between the reactive and non-reactive scat¬ 
tering amplitudes for the H -b H 2 reaction which signif¬ 
icantly alters the angular dependence of the differential 
cross sections (DCSs). 0, @ Unfortunately, state-resolved 
experiments for H -|- H 2 are very difficult in practice and 
Mead’s predictions have not yet been verified. Though 
the isotopic variants are more accessible experimentally, 
theoretical calculations showed negligible GP effects for 
a wide range of collision energies. [91417j| Some relatively 
small rapidly varying oscillations in the DGS due to the 
GP have been seen in the theoretical DGSs at energies 
below that of the GI. At energies above the 

GI, large GP effects on the DGS’s were observed which 
give rise to broader bi-modal features. [Tsl - f^ However, 
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GP effects remained elusive in the integral cross sections 
or reaction rates at any energy. A recent experimental 
attempt to measure the GP oscillations in the DGSs for 
the H -b HD —^ H -b HD reaction at energies below the 
GI was unsuccessful. M 

Until recently, [23 all previous theoretical predictions 
of GP effects on chemical reactivity and experimental 
attempts at its detection have been done at thermal en¬ 
ergies. Recent experimental progress in the cooling and 
trapping of molecules presents a novel energy regime at 
sub-Kelvin temperatures to explore GP effects in chemi¬ 
cal reactions. |23l l24l| In the zero-temperature limit where 
only s-wave contributes, the reaction rates obey the well 
known Bethe-Wigner threshold laws and approach finite 
measurable values. [^.[^ In this Letter, it is shown that 
due to the unique properties associated with ultracold 
collisions, namely: (1) isotropic scattering, and (2) quan¬ 
tized scattering phase shifts, the maximal possible inter¬ 
ference between the different scattering pathways around 
a GI becomes possible. The maximal interference ef¬ 
fects are shown to occur in the fundamental hydrogen 
exchange reaction which results in very large GP effects 
effectively turning on or off the reactivity. The quantized 
scattering phase shifts are a general property of ultracold 
collisions and can occur for interaction potentials which 
support bound states as well as those which do not (as 
demonstrated in this Letter for the H 3 system). In the 
latter case, suitable vibrational excitation of the reactant 
diatomic molecule is required which results in an effective 
reaction pathway (along the vibrational adiabat) that is 
barrierlessfrzj and exhibits a potential well.fH-js^ 

Eigure [T] plots a 2D slice of the 3D ground state H 3 
electronic PEsfall and reaction pathways for the D -b 
HD —>■ D -b HD (panel a) and D -b HD —>■ H -b D 2 (panel 
b) reactions. Hyperspherical coordinates are used which 
have the advantage of showing all arrangement channels 
simultaneously as well as the prominent GI located near 
the center of the plot.js^ Eigure [U corresponds to a stere¬ 
ographic projection of the upper half of the hypersphere 
with a fixed hyperradius of p = 3.75 Uq. The zero of en¬ 
ergy is the bottom of the asymptotic H -b D 2 potential 
well. The contour lines are separated by 2,900K ex- 
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FIG. 1. (color online) A 2D slice of the 3D Born- 

Oppenheimer PES for the HD 2 system is plotted at a fixed 
hyperradius p — 3.75 Uq. The different scattering pathways 
around the Cl are indicated, panel (a) for D + HD —>■ D + 
HD, and panel (b) for D + HD —^ H + D 2 . For clarity, path¬ 
ways for only one of the symmetric HD channels is depicted 
(see text for discussion). 


cept for the two closely spaced contours at 4,640 K and 
5,220K. For clarity, a cut-plane is used at 33, 640 K so 
that the extremely repulsive regions for each channel are 
not plotted. The energy of the Cl in Fig. [T] is 37, 700 K 
and the PES for the excited electronic state is not plot¬ 
ted. 

The two panels in Figure [T] depict the interference 
pathways which can lead to significant GP effects for each 
reaction. In general the total scattering amplitude can 
be decomposed into contributions from each pathway la¬ 
beled by the 1, 2, and 3 in panels (a) and (b). [^. 

For the inelastic scattering in[T](a), pathway #1 (black) 
corresponds to a non-reactive process and pathways ^2 
and #3 (red) corresponds to an exchange process where 
the two identical D nuclei in each HD channel are ex¬ 
changed. For the reactive scattering in [T] (b), pathway 
#1 (black) corresponds to a direct reaction process and 
pathway #2 (red) corresponds to a looping reaction pro¬ 


cess. For the D-atom exchange in [T] (a), it has been 
shown that, due to the direct collinear nature of the re¬ 
action, contributions from pathway ^3 (thin red curve) 
are negligible even for high collision energies approach¬ 
ing that of the CL nn Thus, the total scattering am¬ 
plitude which does not include the GP can be written 
as = (l/V2)(/i + h) where NGP denotes “No 

Geometric Phase” and /i and /2 are the scattering am¬ 
plitudes for pathways #1 and #2 in Fig. [T] (a). [^. [l 3 - [T^ 
The GP alters the sign on the scattering amplitude for 
pathway #2 across the branch cut (black dashed curve 
in Fig. [1]). Thus, the total scattering amplitude which 
includes the GP is given by = (l/\/ ^)( fi -- f?) where 
GP denotes “with Geometric Phase”. [j fl4jl6l | The same 
expressions hold for the NGP and GP scattering ampli¬ 
tudes in Fig. [T](b).[l^-[l^ The encirclement of the Cl by 
the combined pathways #1 and #2 is obvious in Fig. [1] 
(b) but not so obvious in Fig. [I](a). Pathway #2 in Fig. 
[T] (a) encircles the Cl through “symmetric encirclement” 
(i.e., via the symmetrization of the wave function with re¬ 
ject to permutation of the identical D nuclei). sain- 

Q 

The cross sections and rate coefficients are computed 
from the square modulus of the total scattering ampli¬ 
tude ll/ll = (1/2) (// -b /I ± 2 /i /2 cos A) where the -b 
and — denote NGP and GP, respectively. The complex 
scattering amplitudes are expressed as /i = /i exp(f (5i), 
/2 = /2 exp(i < 52 ) and the phase difference A = (52 — <5i. If 
the square modulus of the scattering amplitude for one 
of the pathways is much larger than the other: // >> /| 
or /I >> ff, then the square modulus of the total scat¬ 
tering amplitude is given by ||/|| « ///2 or ||/|| r; /|/2, 
respectively, and the GP effect is negligible. However, 
when the squared moduli are similar // « /|, then 
ll/ll ~ /^(licosA) where / = /i si / 2 . Thus, 
depending upon the sign and magnitude of cos A, the 
reactivity can be dramatically enhanced or suppressed. 
The maximum interference occurs when |cosA| = 1. If 
±cosA = -bl then maximum constructive interference 
occurs and ||/|| ~ 2/^, whereas for ±cosA = — 1, maxi¬ 
mum destructive interference occurs and ||/|| ~ 0. That 
is, if I cos AI ft! 1 then the reactivity can be turned on or 
off by the sign of the interference term. Since the GP al¬ 
ters the sign of the interference term, the GP effectively 
controls the reactivity. 

The quantum reactive scattering calculations for the 
H 3 system were done using a numerically exact time- 
independent coupled-channel method based on hyper- 
spherical coordinates and the GP effect is included us¬ 
ing the vector potential approach, (l^ The com¬ 

puted scattering {S) matrices include all open reactant 
and product diatomic vibrational and rotational states 
on a grid consisting of 112 collision energies spanning 
the range from 1 fiK to 4,640 K relative to asymptotic 
energy of HD(u = 4, j = 0) for the H -b HD and D 
-b HD reactions, and H 2 (u = 4, j = 0) for H -b H 2 . 
For the highly excited reactant vibrational states, the 
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reaction becomes effectively barrierless and exhibits sig¬ 
nificant reactivity at ultracold collision energies, The 
asymptotic energies for HD('(; = 4, j = 0) and H 2 (u = 4, 
j — 0) are 22,109 K and 25, 078 K relative to the bottom 
of the asymptotic diatomic potential wells, respectively. 
The scattering calculations were carried out using two 
accurate ab initio electronic PESs for the H 3 system: the 
BKMP2[3l| and the newer one by Mielke, et al.js^ 



FIG. 2. (color online) The total reaction rate for the H + 
H 2 (i’ = 4, j = 0) ^ H -I- H 2 (para-para) reaction is plotted 
as a function of collision energy (see text for discussion). 
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FIG. 3. (color online) Reaction rates for the D -I- HD(n = 4, 
J = 0) —>■ D -I- HD(a', j') reaction are plotted as a function 
of collision energy: (a) v' = 0, j' — 11, (b) v' = 3, j' = 5, 
(c) v' — 1, and (d) total. The results are for even exchange 
symmetry and include all values of total angular momentum 
J — 0 — 2 (see text for discussion). 

Figure [2] plots the total reaction rate for H -|- H 2 (u = 4, 
j = 0) —1 H -|- H 2 summed over all product vibrational 
and even rotational states (i.e., the para-para transi¬ 


tions). The results which include (do not include) the 
GP are plotted in red (black). The thick curves include 
all values of total angular momentum J = 0 — 2 and 
the thinner curves are for J = 0 only. The rates for 
each value of J are well converged over the entire energy 
range. The total rate is well converged with respect to 
the sum over J = 0 — 2 up to about 2K.[35 The dashed 
and long-short dashed curves compare the results based 
on the BKMP2 and Mielke PES, respectively. Both PESs 
give similar results and predict that the GP enhances the 
ultracold reactivity by a full order of magnitude. Figure 
[3] plots several representative rate coefficients for the D 
-I- HD(u = 4, j = 0) — D -I- HD(z;', j') reaction using 
the BKMP2 PES. The rates are computed for even ex¬ 
change symmetry (i.e., the nuclear motion wave function 
is symmetric with respect to permutation of the iden¬ 
tical D nuclei) and include all values of total angular 
momentum J = 0—2. The rotationally resolved rate 
for v' = 0, j' = II shows that the geometric phase re¬ 
duces the reactivity by over four orders of magnitude. 
The rate for v' = 3, j' = 5 and the vibrationally re¬ 
solved rate for v' = 1 are reduced by nearly two orders 
of magnitude when the GP is included. The total rate 
which includes the GP is reduced by a factor of 50. The 
rates for odd exchange symmetry are similar in magni¬ 
tude except that the GP increases the reactivity in this 
case. [ 3 ^ When both the symmetric and antisymmetric 
nuclear spin states of D 2 are present, the scattering re¬ 
sults must be summed over both even and odd exchange 
symmetries including the appropriate nuclear spin sta¬ 
tistical weights 2/3 and 1/3, respectively. The GP effect 
is reduced but remains significant (see Table U and Ref. 

[si). 

Table U lists a representative sample of the ultracold 
(Ip,K) reaction rates computed for the H 3 system both 
with and without the geometric phase and different ex¬ 
change symmetries. Most notable are the very large GP 
effects seen in the rotationally resolved rates for the D 
-I- HD — >■ D -I- HD reaction for each exchange symmetry. 
These differences persist, albeit smaller when summed 
over both exchange symmetries. The GP effects for the 
H -I- HD — >■ H -I- HD reaction are overall smaller than 
those for the D -|- HD reaction. For the H -|- H 2 —>■ H 
-|- H 2 para-para reaction, the results have already been 
summed over the appropriate nuclear spin states. Large 
differences (~ lOx) between the GP and NGP rates are 
observed even when summed over all v' and j' . The total 
rates summed over all v' and j' using the PES of Mielke 
et al. are tabulated in the last column. They are very 
similar at the state resolved level as well but the over¬ 
all reactivity is slightly reduced (see Fig. [ 21 ). For the 
D -1- HD —H -I- D 2 and H -|- HD —>■ D -|- H 2 reactions 
(not tabulated), the GP and NGP rates are nearly identi¬ 
cal even at the rotationally resolved and single exchange 
symmetry level. The same applies to the H -|- H 2 para- 
ortho reaction (not tabulated). The lack of GP effects 
for these ultracold reactions is due to the direct collinear 
nature of the reaction which results in a tiny contribution 
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TABLE 1. Ultracold (1 /iK) reaction rates for the X + 
HD(?; = 4, j = 0) ^ X + HD(w', f) with X=D and H, 
and the H + H 2 (n = 4, j = 0) —>■ H + H 2 (i>', j' even) re¬ 
actions. The “evn” and “odd” denote exchange symmetry 
and the “GP” and “NGP” denote the results computed with 
and without the geometric phase, respectively. The “(evn) 
-I- (odd)” denote the summed results over even and odd ex¬ 
change symmetry. All rates include the appropriate nuclear 
spin statistical weights and are in cm®/s. 


Reaction 

U=0,i'=ll 

v'=2,Y^., Total 

Mielke PES 

D -f HD (evn) GP 
NGP 

D -b HD(odd) GP 
NGP 

1.01(-18) 

1.13^14) 

5.48^15) 

1.90(-19) 

1.51(-15) 1.05(-14) 
1.52(-13) 4.94^13) 
7.52(-14) 2.44(-13) 
7.42(-16) 4.99^15) 

6.64(-15) 

2.51^13) 

1.25^13) 

2.97^15) 

(evn) -b (odd) GP 
NGP 

5.48^15) 

1.13^14) 

7.67^14) 2.54^13) 

1.53^13) 4.98^13) 

1.32^13) 

2.54^13) 


v'=0,j'=2 

F=0,^,, 


H -b HD (evn) GP 
NGP 

H + HD (odd) GP 
NGP 

1.60(-17) 

1.22(-14) 

3.79^14) 

1.95^17) 

1.58(-14) 9.87(-14) 
8.80^14) 3.15^13) 
2.71^13) 9.35^13) 
4.34(-14) 2.85^13) 


(evn) -b (odd) GP 
NGP 

3.79(-14) 

1.22(-14) 

2.87^13) 1.03^12) 

1.31^13) 6.00^13) 



v'=3,f=4 



H -b Ha GP 

NGP 

4.32(-ll) 

2.16^13) 

8.39(-12) 2.48(-ll) 
5.41(-13) 2.02(-12) 

2.16(-11) 

1.76^12) 


from the scattering amplitude corresponding to looping 


tentials corresponding to the two reaction pathways in 
Fig.m: a) differ by an even (odd) number, then cos A = 1 
(cos A = —1). Maximum constructive (destructive) in¬ 
terference will occur between the two scattering ampli¬ 
tudes contributing to and the opposite interfer¬ 
ence behavior will occur for Thus, the unusually 

large GP effects reported here originate from the isotropic 
(s-wave) scattering and the effective quantization of the 
scattering phase shift which results in | cos A| r; 1.[^I^ 
The mechanism is general and is expected to hold for 
many molecules which exhibit CIs and for which the PES 
and/or the choice of reactant and product states allows 
for a favorable encirclement. 

We emphasize that the interference mechanism re¬ 
ported here is a general property of ultracold collisions 
and will also occur in molecules without CIs or GP ef¬ 
fects. In general, large interference effects can be ex¬ 
pected for barrierless reaction paths which proceed over 
a potential well (due to the PES or vibrational excitation) 
and include contributions from two interfering pathways 
(such as reactive and non-reactive). Experimentally the 
interference (and hence reactivity) might be controlled 
by the selection of a specific nuclear spin state or by the 
application of external electric or magnetic fields to (1) 
alter the relative number of bound states in the effective 
potential wells along each interfering pathway, or (2) alter 
the relative magnitude of the two interfering scattering 
amplitudes. 


pathway in Fig. [T] (b). 


For ultracold collisions of H or D with a high vibra- 
tionally excited HD or H 2 diatomic molecule, leading to 
vibrational quenching, the reaction pathway is effectively 
barrierless with an attractive potential well. [iS-lli Thus, 
each scattering pathway in Fig. [T]can be represented by 
a simple ID spherical well model. For this ID model, 
the scattering phase shifts are known analytically and in 
the zero energy limit they become effectively quantized 
(i.e., they approach mr where n denotes the number of 
bound states in ID spherical well). [l^, If the number 
of bound states in the two different ID spherical well po- 


Acknowledgments 

BKK acknowledges that part of this work was done under 
the auspices of the US Department of Energy under Project 
No. 20140309ER of the Laboratory Directed Research and 
Development Program at Los Alamos National Laboratory. 
Los Alamos National Laboratory is operated by Los Alamos 
National Security, LLC, for the National Security Adminis¬ 
tration of the US Department of Energy under contract DE- 
AC52-06NA25396. The UNLV team acknowledges support 
from the Army Research Office, MURI grant No. W911NF-12- 
1-0476 and the National Science Foundation, grant No. PHY- 


1205838. 


1. H. C. Longuet-Higgins, U. Opik, M. H. L. Pryce, and 
R. A. Sack, Proc. R. Soc. London, Ser. A 244, 1 (1958). 

2. G. Herzberg and H. G. Longuet-Higgins, Discuss. Fara¬ 
day Soc. 35, 77 (1963). 

3. G. A. Mead and D. G. Truhlar, J. Chem. Phys. 70, 2284 
(1979); 78, 6344E (1983). 

4. G. A. Mead, Chem. Phys. 49, 23 (1980). 

5. M. V. Berry, Proc. R. Soc. London Ser. A 392, 45 (1984). 

6. A. Bohm, Quantum Mechanics: Foundations and Appli¬ 
cations, 3rd Edition, (Springer-Verlag, New York), Chap¬ 
ters XII and XIII, (1993), (2001). 

7. G. A. Mead, J. Chem. Phys. 72, 3839 (1980.) 

8. B. Lepetit, A. Kuppermann, Chem. Phys, Lett. 166, 581 
(1990). 

9. B. K. Kendrick, J. Chem. Phys. 112, 5679 (2000); 114, 


4335E (2001). 

10. B. K. Kendrick, J. Chem. Phys. 114, 8796 (2001). 

11. B. K. Kendrick, J. Chem. Phys. 118, 10502 (2003). 

12. B. K. Kendrick, J. Phys. Chem. 107, 6739, (2003). 

13. C. Juanes-Marcos, S. C. Althorpe, J. Chem. Phys. 122, 

204324 (2005). 

14. J. C. Juanes-Marcos, S. C. Althorpe and E. Wrede, Sci¬ 
ence 309, 1227 (2005). 

15. S. C. Althorpe, J. Chem. Phys. 124, 084105 (2006). 

16. S. C. Althorpe, T. Stecher, and F. Bouakline, J. Chem. 
Phys. 129, 214117 (2008). 

17. F. Bouakline, Chem. Phys. 442, 31-40 (2014). 

18. F. Bouakline, B. Lepetit, S. C. Althorpe, and A. Kup¬ 
permann, in “The Jahn-Teller effect fundamentals and 
implications for physics and chemistry”, Castleman, Jr., 












5 


A. P., Toennies, J. P., Yamanouchi, K., & Zinth, W. Eds. 
Springer Series in Chemical Physics, 97, 201-237 (2009). 

19. F. Bouakline, S. C. Althorpe, and D. P. Ruiz, J. Chem. 
Phys. 128, 124322 (2008). 

20. F. Bouakline, S. C. Althorpe, P. Larregaray, and L. Bon¬ 
net, Molec. Phys. 108, 969-980 (2010). 

21. J. Jankunas, M. Sneha, R. N. Zare, F. Bouakline, and S. 
C. Althorpe, J. Chem. Phys. 139 , 144316 (2013). 

22. B. K. Kendrick, J. Hazra and N. Balakrishnan, Nat. 
Commun. in press. 

23. L. D. Carr, D. DeMille, R. V. Krems, J. Ye, New. J. 
Phys. 11, 055049 (2009). 

24. S. Ospelkaus et al., Science 327 , 853 (2010). 

25. H. A. Bethe, Phys. Rev. 47 , 747 (1935). 

26. E. P. Wigner, Phys. Rev. 73, 1002 (1948). 

27. I. Simbotin, S. Ghosal, and R. Cote, Phys. Chem. Chem. 
Phys. 13 , 19148 (2011). 

28. W. H. Miller, N. C. Handy, and J. E. Adams, J. Chem. 
Phys. 72 , 99 (1980). 

29. D. G. Truhlar, B. C. Garrett, and J. S. Klippenstein, J. 
Phys. Chem. 100 , 12771 (1996). 

30. J. Jankunas, M. Sneha, R. N. Zare, F. Bouakline, S. 
C. Althorpe, D. Herraez-Aguilar, and F. J. Aoiz, PNAS 
111 , 15 (2014). 

31. A. I. Boothroyd, W. J. Keogh, P. G. Martin, and M. R. 
Peterson, J. Chem. Phys. 104 , 7139 (1996). 

32. R. T Pack, G. A. Parker, J. Chem. Phys. 87, 3888 (1987). 

33. B. Kendrick and R. T Pack, J. Chem. Phys. 104, 7475 
(1996). 

34. B. K. Kendrick, R. T Pack, R. B. Walker and E. F. Hayes, 
J. Chem. Phys. 110, 6673 (1999). 

35. C. Makrides et al., Phys. Rev. A 91 , 012708 (2015). 

36. See the online supplementary material for more details. 

37. S. L. Mielke, B. C. Garrett, and K. A. Peterson, J. Chem. 
Phys. 116 , 4142 (2002). 



6 


SUPPLEMENTARY MATERIAL 

The three-body quantum reactive scattering methodology is based on a time-independent coupled-channel approach 
using hyperspherical coordinates [3^ and is well suited for treating ultracold collisions, The method is numerically 
exact and accurately treats the body-frame Eckart singularities [3^ associated with non-zero total angular momen¬ 
tum J and includes the geometric phase using the general vector potential approach. [nisi A brief summary of the 
methodology is given here. In the interaction region where the three atoms are in close proximity (i.e., for small 
hyperradius p), Smith-Johnson symmetrized hyperspherical coordinates are used. For larger hyperradius where the 
reactant and product channels become well defined, a properly symmetrized set of Fock-Delves hyperspherical coor¬ 
dinates are used (one for each channel). The hyperradius p is discretized into a set of sectors spanning the range from 
small to large p. The three-body Hamiltonian is diagonalized at each fixed value of p to obtain a set of 2D angular 
wave functions. The 2D angular solutions are independent of the collision energy so only have to be computed once 
for each value of total angular momentum J and inversion parity. The 2D angular solutions form the basis set for 
the coupled-channel equations in p and are used to compute a set of potential coupling matrices within each sector 
and the overlap matrices between the different 2D solutions at the boundaries of each sector. The coupled-channel 
equations are solved for a given collision energy using Johnson’s log-derivative propagator method from small to large 
p. Finally, the asymptotic boundary conditions are applied at large p to compute the scattering S matrix from which 
the cross sections and reaction rates can be computed. 

The individual contributions from each value of total angular momentum J = 0 - 2 to the total reaction rate for the 
D J- HD(u = 4, J = 0) —>■ D J- HD reaction are plotted in Fig. 31] as a function of collision energy. The results with 
(without) the GP are plotted in red (black). The thick red and black curves correspond to the total rate summed over 
all values of J = 0 — 2. The contributions to the rate from J > 0 rapidly decrease with decreasing collision energy due 
to the angular momentum barrier in the entrance channel which increases with J. The rates for each value of J are 
well converged over the entire energy range. Fig. 311 shows that the total rate summed over J = 0 — 2 is converged 
for collision energies up to about 2 K. 
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FIG. SI. The individual contributions to the total reaction rate for D -|- HD(a = 4, j = 0) —^ D -|- HD from each value 
of total angular momentum J — 0 - 2 are plotted (see legend). The results are for even exchange symmetry. The GP and 
NGP reaction rates summed over all values of J = 0 - 2 are shown as thick red and black curves, respectively (see text for 
discussion). 

The reaction rates for the D + HD(u = 4, j = 0) —>■ D J- HD('(;', j') reaction are plotted in Fig. 31|as a function of 
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collision energy for several values of v' and j' (the same as those in Fig. [3] of the main article). The rates correspond 
to odd exchange symmetry and include all values of total angular momentum J = 0 — 2. In contrast to the even 
exchange symmetry results plotted in Fig. |3]of the main article, it is the GP ultracold reaction rates which are largest 
for odd exchange symmetry. Figure 321 plots the reaction rates summed over both even (Fig. [3]) and odd (Fig. 311) 
exchange symmetries multiplied by the appropriate nuclear spin statistical weights. The differences between the GP 
and NGP ultracold reaction rates are smaller but still significant. 





Energy (K) Energy (K) 

FIG. S2. Reaction rates for the D + HD(i; = 4, j = 0) —>■ D + HD(w', j') reaction are plotted as a function of collision 
energy: (a) v' = 0, j' = 11, (b) v' = 3, j’ = 5, (c) v' = 1, and (d) total. The results are for odd exchange symmetry and include 
all values of total angular momentum J = 0 — 2. 


Figure 321panel (a) plots the ratio of the average square moduli of the exchange and non-reactive pathways depicted 
in Fig. [T]of the main article as a function of collision energy for the D -|- HD(n = 4, j = 0) —?> D -f HD('(;' = 3, j' = 5) 
reaction. The results for J = 0, 1, and summed over J = 0 —2 are plotted in black, red, and blue, respectively. In each 
case, the ratio is near unity for a wide range of energies from 1 /iK to 100 K. This implies that maximal interference 
between the exchange and non-reactive pathways depicted in Fig. [T] (a) is possible, and that the interference will be 
governed by the sign and magnitude of cos A. Figure 321 panel (b) plots the average value of cos A as a function of 
collision energy for rriji = 5 (the results for other values of niji are similar). The results for J = 0 and 1 show that the 
scattering phase difference A lies near an even and odd value of tt, respectively. This leads to maximal interference 
between the non-reactive and exchange pathways in Fig. [T](a). Since cos A = -fl for J = 0, the square moduli of the 
total NGP and GP scattering amplitudes are given by ||/'^^^|| = (l-|-cos A) « 2 and ||/^^|| = (I —cos A) ^ 0 
which explains why the NGP ultracold reaction rates are much larger than the GP ones in Fig. [3] of the main article. 
For J = 1 the situation is reversed since cos A = — I. This dramatically reduces the differences between the NGP and 
GP reaction rates at higher collision energies when the contributions from J = 1 and higher J become important. As 
shown in Fig. 321 panel (b), at the higher collision energies (i.e., > 1 mK) where more values of J contribute (blue 
curve), the cos A is no longer quantized and tends to oscillate about zero. For odd exchange symmetry, the situation 
is reversed and cos A = — 1 and -1-1 for J = 0 and 1, respectively. This explains why the GP ultracold reaction rates 
are much larger than the NGP ones in Fig. 33 

The tendency for cos A to be near a multiple of tt (as discussed above in the context of Fig. 321 ) at ultracold 
collision energies occurs for the majority of product states in all three exchange reactions. Figure 33 plots cos A vs 
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FIG. S3. Reaction rates for the D + HD(ii = 4, ji = 0) —>■ D + HD(w', j') reaction are plotted as a function of collision 

energy: (a) v' = 0, j' = 11, (b) v' = 3, j' = 5, (c) v' = 1, and (d) total. The results are summed over both even and odd 
exchange symmetry and include all values of total angular momentum J = 0 — 2. 


COS A at the ultracold collision energy of 1 /iK for the three exchange reactions: (a) D + HD(z; = 4, j = 0) — ?► D + 
HD(n', /), (b) H + HD(n = 4, j = 0) —H + HD(n', j'), and (c) H + H 2 (u = 4, j = 0) —H + H 2 (n', j' even). The 
preference for cos A r; ±1 is striking. The value of cos A for the large majority of product states are clearly clustered 
at the top-right and bottom-left corners. In panels (a) and (b), the values cos A = -1-1 and —1 are preferred by the 
even (red) and odd (black) exchange symmetries, respectively. The quantization effect is most dramatic for the D -|- 
HD reaction (panel a) and the H -|- H 2 reaction (panel c) which explains why these two reactions exhibit the largest 
overall GP effects relative to the H -h HD reaction (panel b) (see also Table U in the main article). In panel (c), the 
cos A are plotted at three different collision energies: 1 (black), IK (blue), and 100/iK (red). At higher collision 
energies the cos A are not quantized and take on a wide range of values with | cos A| < 0.8. However, at the ultracold 
energy of 1 /iK (black), the cos A become essentially quantized and approach values very close to -1-1. 

The dynamical origin of the phase quantization observed in Figs. ^ and ^ can be understood by considering 
a simple ID spherical well potential. Figure 31] shows a ID spherical well potential with depth Vo and radius To- 
The scattering phase shift S can be derived analytically for this model and is given by the expression shown in Fig. 
m for s-wave scattering (i.e., 1 = 0). The analytic phase shift is plotted as a function of the wave number k and 
each blue curve corresponds to a different well depth Vo- In the ultracold limit fc —^ 0, the majority of the phase 
shifts (blue curves) converge towards a multiple of tt. For very small values of k, the phase shift quickly jumps from 
one value of tt to another as the well depth Vo is varied. The jumps occur at intervals of tt/2 which correspond to 
resonances associated with a continuum state dropping into the well to become a bound state. In general, the phase 
shift for fc —>■ 0 is equal to n tt where n denotes the number of bound states supported by the potential well. For 
high-lying vibrationally excited reactant states, the ultracold hydrogen exchange reaction is effectively barrier less ( 2 ^ 
and proceeds along an effective attractive potential well (i.e., vibrational adiabatl. f^sl - lsoll Thus, the scattering phase 
shifts along each pathway depicted in Fig. [T]of the main article can be modeled using an effective ID spherical well 
potential. Since each pathway has a different effective well depth Vo and radius Tq, it will support a different number 
of bound states and hence the ultracold scattering phase shift along each pathway will approach a different multiple 
of TT. This implies that the phase difference A will also approach a multiple of tt. If the multiple of tt is an even or odd 
number, then cos A = -|-1 or cos A = —1, respectively. This ultimately determines whether the interference will be 
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FIG. S4. In panel (a), the ratio of the average square modulus of the exchange (pathway ^2 in Fig. [l]a) and non-reactive 
(pathway in Fig. [T]a) contributions to the total scattering amplitnde is plotted as a function of collision energy for the D 
-I- HD(u = 4, 7 = 0) D -f- HD(u' = 3, j' = 5) reaction for even exchange symmetry. The black and red curves correspond to 
single values of J = 0 and 1, respectively. The blue curve corresponds to the sum over all values of J between 0 and 2. Panel 
(b) plots the average of cos A as a function of collision energy for m^/ = 5 using the same color designations as in panel (a). 
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FIG. S5. The cos A vs cos A are plotted for three hydrogen exchange reactions. In panel (a), the results for D -|- HD(u = 4, 
7 = 0) —>■ D -b HD(u', j') are plotted at 1 fiK. Each point corresponds to a v', j' product state and red (black) points correspond 
to even (odd) exchange symmetry. Panel (b) is similar to panel (a) except that the results are for H -|- HD(w = 4, 7 = 0) —>■ H 
+ HD(u', j'). Panel (c) plots the results for H + H 2 (u = 4, 7 = 0) ^ H -|- H 2 (^)^ j' even) at three different collision energies: 
black (1/rK), blue (IK), and red (100 K). 


constructive or destructive and whether the GP or NGP reaction rates will be the largest. Furthermore, the effective 
ID spherical well potential associated with each of the reaction pathways also depends upon the particular reactant 
and product states involved (hence the distribution of points observed in Fig. ^S]). 

































10 


V(r) 

0 


Spherical Well Potential 


Scattering Phase Shift 

k cos(kro)sin(k'rJ - k'cos(k'ro)sin(kr„) 
k'cos(k'ro)cos(krJ + k sin(k'rjsin(kr„) 


5/71 



FIG. S6. The ID spherical well potential model with depth K> and radius ro, and its associated scattering phase shift 
(for I — 0) are plotted. The scattering phase shifts 5 in units of tt are plotted as a function of the wave number in units 
Tv/ro for different potential well depths Vo = q^/{2fi) separated hy q — 0.025 rr/ro (blue curves) where k'^ — q^ = . In the 

zero temperature limit fc —^ 0, the majority of scattering phase shifts approach a multiple of tt (i.e., they become effectively 
quantized). 































